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FOREWORD 


This  report  is  for  Phase  IV  of  the  program  entitled  "The  Effect  of  Envi¬ 
ronment  on  Mechanical  Behavior  of  AS/3501-6  Graphite /Epoxy  Material."  The 
work  was  conducted  at  Vought  Corporation  Advanced  Technology  Center.  This 
report  covers  all  contract  work  in  the  period  3  August  1981  to  3  November  1982 
and  was  sponsored  by  the  Naval  Air  Systems  Command  under  Contract  Number 
N0001S-81  -C-0179. 

Mr.  M.  Stander  was  the  Navy  Project  Manager  and  Mr.  R.  C.  Knight  was  the 
Vought  Project  Manager.  Key  technical  personnel  were  Dr.  T.  Ho,  Principal  In¬ 
vestigator,  and  Dr.  R.  A.  Schapery ,  Texas  A&M  University,  Technical  Consultant 
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1.0  INTRODUCTION 


Wide  usage  of  advanced  composite  materials  in  future  aerospace  systems  is 
projected  to  become  a  reality  over  the  next  few  years.  A  V/STOL  fighter  air¬ 
craft  is  one  of  several  Navy  articles  expected  to  be  increasingly  dependent  on 
composite  materials  for  primary  and  secondary  structural  components. 

Presently,  the  response  of  advanced  composite  materials  to  arbitrary  load 
and  environments  is  known  to  have  a  significant  adverse  effect  on  the  materi¬ 
als  performance.  The  materials  reponse  in  severe  environments  and  related 
fatigue  lifetime  prediction  methods  are  presently  insufficient.  It  is 
believed  this  composites  research  program  assists  in  resolving  this  deficiency 
in  the  technology  base. 

The  overall  objectives  of  this  research  program  are: 

o  To  ascertain  if  the  mechanical  response  of  AS/3501-6  graphite/epoxy 
composite  material,  subject  to  various  time,  temperature  and  moisture 
effects,  car.  be  characterized  using  traditional  viscoelastic  shift  fac¬ 
tors,  and  to  formulate  a  master  curve  of  material  property  dependence 
on  time,  temperature  and  humidity  if  possible. 

o  To  ascertain  the  feasibility  of  predicting  fatigue  failure  of  a  compos¬ 
ite  material  by  accounting  for  the  linear  viscoelastic  behavior  of  the 
resin  in  various  temperature  and  humidity  environments. 

o  To  determine  if  a  specific  thermal  conditioning  environment  can  be 
directly  substituted  for  a  specific  moisture  conditioning  environment, 
over  a  prescribed  temperature  vs.  humidity  range  for  AS/3501-6  graph¬ 
ite/epoxy  material,  and  obtain  an  equivalent  moisture  effect  on  mechan¬ 
ical  and  fatigue  properties.  If  this  can  be  shown,  a  substantial  cost 
and  time  savings  in  moisture  conditioning  cf  the  test  specimen  can  be 
achieved  and  possibly  extended  to  other  composite  specimens. 
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Through  the  research  of  the  Phase  I  program1  the  first  and  third  objec¬ 
tives  previously  mentioned  were  accomplished.  Thus,  linear  viscoelast{city  is 
an  appropriate  method  to  characterize  the  composite's  response  with  respect  to 
time,  temperature  and  humidity.  A  definitive  equivalence  relationship  between 
one  environment  and  another  is  not  feasible  in  the  sense  that  each  tempera¬ 
ture/humidity  condition's  effect  on  the  materials  mechanical  response  is 
unique  i n  itsel f . 

The  research  objective  of  this  Phase  IV  program  was  to  continue  the  Phase 
II  and  Phase  III  effort  to  accomplish  the  second  objective,  i.e.,  to  predict 
the  fatigue  failure  of  composite  material  by  accounting  for  its  viscoelastic 
behavior  in  various  temperature  and  humidity  environments. 

The  basic  Phase  IV  program  scope  was  increased  to  evaluate  the  edge  crack¬ 
ing  and  delamination  of  the  (+45/90o)c  layup  to  supplement  the  fatigue 

c  5  p  3 

data  obtained  from  the  previous  Phases  II  and  III.  This  added  scope 
entailed  the  residual  strength  tests  and  delamination  tests  on  graphite/epoxy 
coupons.  Since  the  fatigue  analysis  requires  that  the  state  of  stress  of  the 
lamina  be  known,  off-axis  specimen  tests  were  conducted. 

The  Phase  IV  program  schedule  is  shown  in  Figure  1. 


2 


TASKS 


1.0  SPECIMENS 

1.1  SPECIMEN  FABRICATION 

1.2  SPECIMFN  PREPARATION 

AND  CONDITIONING 

2.0  TESTS 

2.1  FATIGUE  TEST  OF  (+45/ 
902)s  -  SPECIMENS 

2.2  STATIC  TEST  OF  OFF-AXIS 
SPECIMENS 

2.3  FATIGUE  TEST  OF  (90 J 
+45 )  SPECIMENS  £ 

3. U  ANALYSIS 

3.1  DATA  REDUCTION 

3.2  TEST/ANALYSIS  CORRELA¬ 
TION 

3.3  THEORETICAL  DEVELOP¬ 
MENT  FOR  THE  MATHEMATI¬ 
CAL  MODEL 

4.0  REPORTING 

4.1  UUARTERLY  REPORTS 

4.2  FINAL  REPORT 


2.0  MICROMECHANICS 


Results  from  Phase  II  and  Phase  III  programs  indicated  that,  for  the 
matrix  dominated  laminate  such  as  the  (+45o/902)$  layup,  composite 
fatigue  characteristics  can  be  interpreted  better  by  relating  to  the  matrix 
through  micromechanics  as  illustrated  in  Figure  2.  In  this  approach,  laminate 
characteristics  are  first  related  to  the  lamina  properties  which,  in  turn,  are 
related  to  the  matrix  properties.  Since  lamina  properties  are  needed  to 
relate  the  laminate  to  the  matrix,  quadratic  tensor  polynomial  criteria  is 
adapted  to  quantify  the  general  stress  state  of  the  laminate. 

In  fatigue  analysis,  the  edge  delamination  type  failure  mode  is  believed 
to  be  involved  in  the  fatigue  failure.  The  separation  of  delamination  failure 
mode  data  from  normal  matrix  damage  degradation  data  is  needed  before  the 
analysis  procedure  in  Figure  2  can  be  used  properly. 
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The  Tsai -Hahn  micromechanics  relationships  have  been  incorporated  into 
the  LAMINATE  ANALYSIS  computer  program.  The  formulas  for  the  micromechanic 
relationships  are  shown  in  Table  1.  These  relationships  originated  from 
Hashin's  composite  cylinder  assemblage  (CCA)  model.5  Especially,  the  value 
of  transverse  shear  modulus,  Gy,  and  transverse  Young's  modulus,  Ey,  are  upper 
bounds  of  the  CCA  model .  Nevertheless,  the  Tsai-Hahn  relationships  are  simple 
enough  for  engineering  applications. 

The  use  of  the  Tsai-Hahn  relationships  for  advanced  composites  is  not 
strai ghcforward  due  to  the  unknown  nature  of  both  its  fiber  and  matrix  proper¬ 
ties,  especially  in  their  in-situ  state.  On  the  basis  of  measured  modulus 
values  of  unidirectional  qraphite/epoxy  composite,  five  anisotropic  fiber- 
properties  (E^,  E£f,  G-j ^ f »  (’23f>  v  12f^  anrl  two  1S0"tr°Pic  matrix 
properties  (E^,  v  m)  can  be  recovered  without  direit  measurement  of  each 
constituent's  properties,  which  are  difficult  to  obtain  with  current  test 
techniques.  The  procedure  for  finding  these  seven  material  properties  is  a 
combination  of  matrix  property  assumptions  and  parametric  sensitivity  rela¬ 
tionships.  Although  this  is  a  cumbersome  approach  to  obtain  the  material 
properties,  it  is  an  easier  task  than  trying  the  analytical  conversion  of  the 
relationships  in  Table  1  and  the  corresponding  expensive  test  program. 


FIGURE  2.  THE  RELATIONSHIP  OF  LAMINATE  FATIGUE  COMPLIANCE 
TO  MATRIX  FATIGUE  COMPLIANCE 


TABLE  1  FORMULAS  FOR  COMPOSITE  MODULI 


ENGINEERING 

CONSTANT 


P  =  1  (VcPr+  nV  P  ) 
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bulk  modulus 
numerical  constant(s) 
fiber  direction 
transverse  direction 

A  sensitivity  study  was  conducted  to  specify  the  different  combinations  of 
the  seven  material  properties  most  sensitive  to  the  known  unidirectional 
data-  Trie  unidirectional  stiffness  (Q.j)  properties  of  a  lamina  and/or  gen¬ 
eral  laminate  can  be  represented  by  the  parametric  forms  shown  in  Table  2 
wnere  fiber  volume  V^-  is  included  for  analysis.  The  Poisson's  ratio  of  the 
matrix  material  is  assumed  to  be  0.35  in  all  of  these  studies.  The  sensitivi¬ 
ty  relations  Of  Qjj  are  obtained 
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oy  taking  tno  differentials  in  the  function.  Equation  (1)  can  then  be 
used  to  evaluate  the  percent  cnange  of  due  to  percent  change  of  the 
parameters  i n  Tabl e  2. 

A1  ternately,  by  assuming  E[n  =  0.6o  msi,  functional  dependence  of  the 
unidirectional  properties  on  each  of  their  dominant  parameters  are  obtained 
using  the  L/V1INATE  ANALYSIS  computer  routine.  Results  are  shown  in  figures  3 
ana  4.  Axiai  modulus  E-|  ^  of  different  graph  ite/  epoxy  layups  are  obtained 
from  test  results  in  References  1,  2,  and  6  and  are  shown  in  Table  3  for  boln 
room  temperature  environment  and  a  170°F/50^  R.H.  environment.  With  tne  in¬ 
formation  from  figures  3  and  4  and  Table  3,  the  fiber  and  matrix  properties 
were  derived,  and  are  shown  in  Table  4,  from  Equation  (1)  and  Table  4,  the 
sensitivity  matrix  for  unidirectional  lamina  and  (jt/(5/yCT>)  s 
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TABLE  2  PARAMETRIC  DEPENDENCE  OF  MATERIAL  PROPERTIES 


TABLE,  i  LONGITUDINAL  MODULUS  PROPERTIES  FOR  AS/3501-6 


ENVIRONMENTS 


75°F/50%  R.H. 


170°F/50%  R.H. 


SPECIMEN  LAYUP 


RESULTS  FROM 
TSAI -HAHN  RELATIONS 
(MSI  ) 


EXPERIMENTAL  RESULTS 
(MS!) 


18.95 

27S0 


TABLE  4  MATERIAL  DATA 


CONSTITUENTS’  DATA  FOR  Gr/Ep 


Elf  = 

28.6  MSI 

E2f  = 

4.2  MSI 

623f 

=  0.89  MSI 

G1 2f 

=  2.50  MSI 

v12f 

-  0.32 

-  0.65 

E  =  0,66  MSI 
m 

vm  =  0o35 

rnMnncTTr  nnrA  rronM  M T rpnMCrMAM TPQ 

E-j  =  18.82  MSI 
E2  =  1.54  MSI 
G12  =  0.805  MSI 
G23  =  0.512  MSI 
v-j  2  =  0.33 
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laminates  as  shown  in  Table  5  may  be  calculated.  The  significant  information 
obtained  from  the  sensitivity  matrix  assuming  that  the  matrix  properties  are 
fixed  is: 

1.  Fiber  volume  fraction  is  the  most  dominant  parameter  affecting  material 
property  change  of  either  a  lamina  or  laminate  (Table  5). 

2.  Unidirectional  properties  E-^,  E^^.  612*  r,23  ancl  v  12  dre 
also  strongly  influenced  by  fiber  properties  Elf,  E 2f,  G-^p, 

and  v -) 2f  respectively. 

These  results  are  helpful  in  adjusting  values  in  Table  4  to  meet  different 
graphite/epoxy  test  results.  With  the  LAMINATE  ANALYSIS  computer  program,  the 
effect  of  matrix  modulus  on  laminate  Ej-|  property  can  be  evaluated.  The 
result  is  shown  in  Figure  5  where  the  compliance  curves,  instead  of  modulus, 
are  compiled.  These  curves  are  important  in  analyzing  both  fatigue  strength 
and  fatigue  failure  data  according  to  the  compliance  relationship  in  Figure 
2.  There  is  one  question  in  correlating  the  LAMINATE  ANALYSIS  results  to  the 
real  test  data.  In  comparing  the  laminate  data  in  Table  3,  especially  those 
from  complicated  lay-ups,  values  of  laminate  E-j  ’  obtained  from  th,(  LAMINATF 
ANALYSIS  program  are  higher  +han  those  obtained  from  tests.  This  situation 
has  been  analyzed  by  Arenhurn7  who  attributes  at  least  some  of  the 
difference  to  transverse  shear  property  effects  as  defined  in  his  finite 
element  analysis.  There  are  more  discussions  on  this  in  Section  4.0. 
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TABLE  5 


SENSITIVITY  MATRIX  FOR  UN  1 0 1  RECT I CNAL  LAM  It* 


AE.f 

AE2f 

AG23f 

mm 

E2f 

G23f 

AE11 

0.988 

o 

o 

En 

^22 

E?2 

0 

0.292 

-0.02& 

AG23 

o 

o 

0.43 

G23 

AG  .2 

G12 

0 

0 

0 

Av12 
"l  2 

0 

0 

0 

SENSITIVITY  MATRIX  FOR  (+4 5/902)s  LAMINATE 


AE,f 

AE2f 

!  ” 

ou23f 

E1f 

E2F 

m 

AEn 

Ei  i 

0.554 

0.083 

-0.012 

AE22 

E22 

0.355 

0.039 

0.003 

“23 

G23 

- 

- 

- 

AG,  2 
"G12 

0.824 

0.011 

0 

Av,2 

"12 

0.050 

0.050 

.0.050 

MATRICES 
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LAMINATE  COMPLIANCE 


3.0  SPECIMEN  FABRICATION  AND  CONDITIONING 


After  reviewing  the  fatigue  data  of  the  AS/3501-6  graphite/epoxy  speci- 
2  3 

mens,  ’  it  seemed  important  to  study  the  influence  of  matrix  degradation  on 
the  edge  del  ami  nation  failure  mode  in  the  fatigue  process.  The  test  program 
was  defined  as  shown  in  Table  6.  The  residual  strength  tests  and  delamination 
tests  are  designed  to  evaluate  the  edge  cracking  phenomenon. 

Hercules  AS1/3501  -6  graphite/epoxy  prepreg  tape  was  used  to  prepare  the 
specimens.  A  total  of  232  specimens  to  be  tested  according  to  test  matrix, 
Table  6,  were  obtained  from  fourteen  panels  of  different  sizes  as  shown  in 
Table  7  and  Figure  6.  The  various  sizes  of  the  panels  resulted  from  the  ply 
lay-up  variation,  defect  inserti  ns  and  the  number  of  specimens  required. 
Typical  specimen  dimensions  are  7"  long  by  0.75"  wide  as  shown  in  Figure  7. 

The  delamination  specimens  have  brass  shim  inserts  (0.001"  thick)  at  the  edges 
as  shown  in  Figure  8  to  simulate  the  edge  cracks.  Typical  layup  procedures 
and  cure  cycles  are  shown  in  Table  8. 

For  moisture  absorption,  steel  chambers  enclosed  in  a  forced  air  Blue-K 
oven  were  used  to  environmentally  condition  the  specimens  to  the  desired  mois¬ 
ture  levels.  A  saturated  aqueous  solution  of  sodium  bromide  (NaBr)  in  contact 
with  a  solid  phase  of  the  salt  at  170°F  was  used  to  generate  the  50%  rela¬ 
tive  humidity  level.  Sodium  fluoride  (NaF)  was  used  to  generate  the  95%  rela¬ 
tive  humidity  level.  Approximately,  fifteen  days  of  conditioning  was  required 
for  the  eight-ply  specimens  and  eighty-five  days  for  the  twenty-ply  speci¬ 
mens.  The  moisture  distribution  inside  the  specimen  should  reach  a  reasonably 
uniform  state  after  the  above  mentioned  conditioning  period.  All  specimen  end 
tabs  were  protected  with  aluminum  backing  tape  wrapping  and  EA934  coating  to 
prevent  moisture  ingress  during  the  environmental  conditioning  process. 
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TABLE  7  COMPOSITE  PANELS  FOR  THE  PROGRAM 


SPECIMEN 

NO.  OF  SPECIMENS 

SPECIAL  FEATURE 

PANEL 

LAY-UP 

FABRICATED 

OF  THE  SPECIMEN 

A 

(15)20 

18 

(30)20 

18 

(60)20 

10 

B 

(O)o 

6 

0 

C 

(45)20 

18 

D 

(90)20 

10 

E 

(+45/90  ) 

—  z  s 

62 

F 

(902/+45)s 

42 

G-l 

(+45/902)s 

4 

WITH  0.02"  INSERTS  IN 

G-2 

(+'+5/902)s 

4 

THE  MIDDLE  LEVEL 

G-3 

(+45/902)s 

4 

WITH  0.04"  INSERTS  IN 

G-4 

(+45/902)s 

4 

THF.  MIDDLE  LEVEL 

G-5 

G-6 

(+45/30,) s 
(+45/30,) 

—  t-  s 

4 

4 

WITH  0.02"  INSERTS  IN 

THE  QUARTER  THICKNESS 
LEVEL 

G-7 

(+45/90,).. 

4 

WITH  0.04"  INSERTS  IN 

G-8 

—  Zb 

(+45/90 ,) 

—  z  s 

4 

THE  QUARTER  THICKNESS 
LEVEL 

P  PROCEDURE  AMO  CURE  CYCLE  FOR  AS/3501 - 


slowly  to  below  150°F  (Cool  no  faster  than  5°F  per  minute 
down  should  take  approximately  45  minutes) 


4.0  STATIC  TESTS  AND  ANALYSIS 


4.1  STATIC  TESTS 

Off-axis  static  tests  were  conducted  according  to  the  test  matrix  of  Table  b. 
The  six  types  of  unidirectional  specimens  for  off-axis  tests  are  (0°)g, 

(15°)20,  (30o)20,  (45o)20,  (60°) 2q,  and  (90o)20.  Some  static 

tests  were  also  conducted  on  (+45  /90°  )  and  {9(P^/H5°  )  laminates  and  the 

data  used  to  verify  the  quadratic  tensor  failure  criterion. 

Tests  were  performed  using  a  Shorewestern  environmental  test  machine  as 
illustrated  in  Figure  9.  Six  temperature  and  humidity  test  environments 
were:  75°F/dry,  75°F/50%  RH,  75°F/95%  RH,  132°F/50%  RH,  132°F/95% 

RH,  and  170°F/50%  RH.  Extensometers  were  used  to  measure  the  strain  at  the 
two-inch  gage  section  of  the  test  specimen.  For  selected  75°F/dry  environ¬ 
ment  tests,  strain  gages  were  also  used  in  conjunction  with  the  extensometers 
in  order  to  verify  the  calibration  of  the  extensometer  system. 


4.2  ANALYSIS 


Test  data  are  shown  in  Figure  10  for  strength  and  Figure  11  for  modulus. 
The  best  fit  of  lower  and  upper  bounds  for  the  test  data  are  also  drawn  in  the 
two  figures  together  with  the  corresponding  classical  laminate  analysis  (CLA) 
comparison  data  to  be  explained  below.  The  Poisson's  ratio  is  0.32  based  on 
the  (0)g  specimen  test  at  75°F/dry  environment  with  strain  gages.  Tabu¬ 
lated  values  of  the  upper  and  lower  bounds  in  Figures  10  and  11  for  various 
angles  are  shown  in  Table  9.  The  shear  modulu  s  G12  can  be  obtained  by  the 
following  relationship: 


2v’2  ♦  i  ,  i 

F  +  5  5  I 

11  sin  ecos  0  E-^ 


4  .  4 

cos  e  sin  e  \ 

r -  -  r -  ) 

L11  b22 


(2) 


where  E^  is  the  longitudinal  modulus,  E22  is  the  transverse  modulus  and 
E-j' -]  is  the  axial  modulus  of  the  off-axis  specimen.  Each  off-axis  specimen 
has  a  calculated  G^-value  and  the  average  is  0.752  msi  as  derived 
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FIGURE  9.  STATIC  TEST  ARRANGEMENT  AND  DATA  ACQUISITION  SYSTEM 


□  75/50 
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FIGURE  10  STRENGTH  OF  STATIC  SPECIMENS 


□  75/50 


TABLE  9  UNIDIRECTIONAL  SPECIMEN  STATIC  TEST  DATA  SUMMARY 


STRENGTH 

(msi) 

«»8 

(15)20 

(3°)20 

<45>20 

<“>20 

O 

C\J 

O 

<7> 

i 

UPPER  BOUND 

0.216 

0.0415 

0.0224 

0.0164 

C.0130 

0.010 

LOWER  BOUND 

0.205 

0.0315 

0.0151 

0.0104 

0.0078 

.  . 

C.  006 

MODULUS 

(msi) 

(0)g 

(15)2q 

Bi 

(S0)20 

1 

<*»20 

UPPER  BOUND 

20.25 

8.6 

r  1 

3.1 

B 

1.5 

1.4 

LOWER  BOUND 

20.25 

7.1 

lJI. 

B 

1.3 

1.2 

i 
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from  the  upper  bound  test  data  and  0.62  msi  from  the  lower  bound.  From  the 
calculated  average  value,  and  Equation  (2),  the  theoretical  modulus  of 
the  off-axis  angle  specimen  is  obtained.  These  theoretical  modules  are  shown 
in  Figure  11.  The  test  value  and  theoretical  value  show  a  favorable  compari¬ 
son.  The  theoretical  values  a^e  used  later  to  interpret  the  fracture  energy 
and  fatigue  data. 

The  strength  data  in  Figure  10  were  used  to  evaluate  the  parameters  of  the 
Tsai-Wu  quadratic  failure  criterion: 


2 

0 

x 

u 

X 


+  F  13 

xv 

+  F  o 

y  y 


o 

x  y 

=  l 


+  F 


ss 


(3) 


where  x  is  the  longitudinal  direction  and  y  is  the  transverse  direction. 

Since  no  compression  or  bi-axial  tests  were  conducted,  compression  and  shear 
data  from  Tsai-Hahn^  was  used.  By  using  the  upper  bound  values  and  lower 
bound  values  from  Figure  10,  parameters  in  Equation  (3)  can  be  derived.  Cal¬ 
culations  and  results  are  shown  in  Table  10.  The  theoretical  strength  diagram 
from  the  Tsai-Wu  criterion  is  shown  in  Figure  12.  The  theoretical  values 
(CLA)  are  also  shown  in  Figure  10  to  compare  with  test  data.  The  static 
strength  values  of  the  ( +45/902 )s  fall  within  the  theoretical  range  for 
(j45/90.?)s  specimens.  The  Tsai-Wu  criterion  can  be  reasonably  used  for 
fatigue  failure  analysis  if  the  residual  strength  is  not  affected  after 
fatigue  cycling. 

To  study  the  failed  specimen,  a  radiographic  technique  was  developed  for 
inspection  of  graphite/epoxy  specimens.  The  delaminations  present  in  the 
specimens  could  not  be  detected  initially  by  conventional  X-ray  methods.  A 
radiographic-opaque  solution  (Znl^)  was  developed  and  applied  to  the  edges 
of  each  of  the  failed  specimens  with  cotton  swabs.  After  a  fifteen  minute 
"dwell  time"  had  elapsed,  the  specimens  were  re-exposed  with  the  x-ray  instru¬ 
ment.  settings  remaininei  the  same.  The  radiograph  from  the  second  exposure 
(with  solution  applied)  enabled  complete  detection  of  the  delaminated  areas  of 
each  specimen  with  excellent  detail  as  shown  in  Figure  13.  Massive  delamina¬ 
tions  in  the  outside  90°  layers  for  (902/j45)s  type  specimens  were 
observed  all  over  the  length  of  the  specimen.  Del  ami  nations  of  90°  layers 
in  the  (+45/902 )$  are  more  confined  to  the  edges. 
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TABLE  10  FAILURE  CRITERION  PARAMETERS 


STRENGTH 

(msl) 

X 

■ 

X 

Y 

1 

Y 

s 

UPPER  BOUND 

0.216 

0.216 

0.010 

0.030 

0.0135 

LOWER  BOUND 

0.205 

0.205 

0.006 

0.018 

0.0081 

..  . 

TSAI-WU  Q 

JADRATIC  FAILU 

RE  CRITERION 

— 

F 

F 

F 

F 

F 

C 

xx 

xy 

yy 

ss 

X 

y 

(msi-2) 

(msi-2) 

(msi-2 ) 

(ms i -2 ) 

(msi-' ) 

(msi-1 ) 

21.43 

-133.6 

3333 

5489 

0 

66.7 

23.80 

-234.7 

9259 

_ 

15242 

0 

111.1 

X  :  LONGITUDINAL  TENSILE  STRENGTH 

X'  :  LONGITUDINAL  COMPRESSIVE  STRENGTH 

Y  ;  TRANSVERSE  TENSILE  STRENGTH 

Y'  :  TRANSVERSE  COMPRESSIVE  STRENGTH 

S  :  SHEAR  STRENGTH 

F  =  1  c  _  1  p  _  _]  /o  /  p  p 

XX  "  xx'  *  yy  ~  ‘yyTt  xy  ~  "  '  V  xx  yy 


ss 


1 

7 
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FIGURE  13.  A  RADIOGRAPHIC  TECHNIQUE  FOR  GR/EP  SPECIMENS 


5.0  FATIGUE  TESTS  AND  ANALYSIS 


In  the  previous  Phase  II  and  Phase  III  programs,  the  mean  fatigue  strain 
data  and  strain  amplitude  for  the  { +45/902 )s  layup  were  developed  as  shown 
in  Figures  14  through  17.  In  this  Phase  IV  program,  limited  numbers  of 
fatigue  tests  on  (+45/902 )$  specimens  were  conducted  to  provide  baseline 
data  in  comparing  the  original  data  (Phase  II  and  Phase  III  data)  to  current 
data.  The  current  fatigue  tests  of  (902/^-5 ) s  were  conducted  at 
75°F/50%  RH  environments  and  132°F/9b%  RH  to  invetigate  the  effect  of  lay¬ 
up  sequence.  The  S-N  data  is  shown  in  Figure  18  and  strain  histories  are 
shown  in  Figures  19  through  22. 

In  order  to  have  a  fair  assessment  of  the  fatigue  data,  basic  static  test 
data  from  both  the  old  batch  of  prepreg  and  current  batch  of  prepreg  is  sum¬ 
marized  in  Table  11.  Due  to  better  static  strength,  the  55%  fatigue  stress 
level  used  for  the  Phase  IV  program  is  12%  higher  than  that  in  the  Phase  III 
program.  The  same  is  also  true  for  the  50%  stress  level  and  60%  stress 
level.  This  difference  in  the  basic  data  raises  the  question  as  to  which  are 
the  reliable  composite  properties  that  can  be  used  for  micromechanics  or  other 
analysis.  To  understand  the  parameters  that  cause  the  variation  of  the  com¬ 
posite  property,  a  compilation  of  test  data  for  AS-graphi te/epoxy  materials 
was  made  and  the  property  list  is  shown  in  Table  12.  The  disparity  among 
basic  properties  from  different  research  sources  for  the  same  material  prompts 
the  necessity  in  looking  into  the  micromechanics  aspect.  Relationships 
between  fiber,  matrix  and  composite  as  illustrated  in  Figures  3  and  4  can  be 
utilized  to  estimate  the  properties  of  the  fiber  from  the  measured  composite 

data.  The  results  are  shown  in  Table  13  together  with  the  data  obtained  by 
7  8 

Arenburg  and  Crane  and  Adams.  In  the  calculations  for  Table  13, 

Poisson's  ratio  v  23  was  obtained  from  E2  and  G23  in  the  isotropic 
plane.  Based  on  the  results  in  Table  13  and  ini cromechanical  relationships  in 
Figures  3  and  4,  variations  of  composite  test  data  seem  to  come  from  three 
possible  areas:  1)  variation  of  fiber  property  from  batch  to  batch;  2)  varia¬ 
tion  of  matrix  characteristics  from  batch  to  batch;  and  3)  inconsistency  in 
precise  layup  and  cure  processing  which  will  result  in  changes,  from  panel  to 
panel,  of  fiber 
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DIFFERENT  SYMBOLS  FOR  DIFFERENT  SPECIMENS 


MAXIMUM  FATIGUE  STRESS  =  60%  UTS  *  DIFFERENT  SYMBOLS  FOR 

DIFFERENT  SPECIMENS 


CYCLES 


75°F/5G%  RH  for  ( ±  45/902)s 
75°F/50%  RH  for  (902/ ±  45)s 
l32°F/50%  RH  for  (902/±45) 
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FIGURE  18  S-N  DATA  OF  (+45/ 902 ) s  AND  (9C2/+45)s  SPECIMEN 
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FATIGUE  STRAIN  CURVE  OF  (<KL/+45}_  SPECIMEN  AT  75°F/503  RH 
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FIGURE  2C  MEAN  FATIGUE  STRAIN  CURVE  OF  (SC9/+4S)c  SPECIMEN  AT  132°F/95%  RH 


MAXIMUM  FATIGUE  STRESS  =  50%  UTS 


100  1,000  1 0,000  100,000 
NUMBER  OF  CYCLES,  N 

FIGURE  21  STRAIN  AMPLITUDE  OF  (90.,/+45)c  SPECIMEN  AT  75°F/50%  RH 


MAXIMUM  FATIGUE  STRESS  =  50%  UTS 


CTi 

CNJ 

ro 

LlJ 

r— 

lO 

t 

LO 

1 

LO 

1 

o 

CVJ 

1 

t 

T 

LO 

UT) 

LO 

□ 

LO 

uO 

LO 

O 

o 

s 

h- 

<c 

s 

S 

o 

cr> 

cn 

cr> 

u. 

cr> 

CT> 

cr» 

> 
H— 1 

> 

4—4 

:> 
i — i 

2- 

>■ 

* — i 

►— t 

»— * 

LjJ 

UJ 

UJ 

S 

<4 

LU 

UJ 

UJ 

CD 

cd 

CD 

s 

CD 

CD 

CD 

CO 

CO 

CO 

1 — 1 

CO 

CO 

00 

O  <1  □ 


O  <J  □ 


(ni/ni  rt)  3amndw\/  nivhis 


LO  LO 
*Cf 

o  o 
cr>  cr» 


UJ  LlJ 
CD  CD 
CO  00 

O  <j 


41 


FIGURE  22  STRAIN  AMPLITUDE  OF  (90,/+45)„ 
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TABLE  13  MICROMECHANICAL  PROPERTIES  OF  AS-TYPE  FIBER  COMPOSITES 
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volume  content,  volatile  content,  and  void  content.  Based  on  these  observa¬ 
tions,  fatigue  data  and  fatigue  analysis  needs  to  be  associated  with  the 
individual  batch  properties.  For  present  analysis  purposes,  the  data  in  the 
last  column  of  Table  13  were  used. 


6.0  RESIDUAL  STRENGTH  TESTS  AND  DELAMINATION  TESTS 


Residual  strength  tests  and  del  ami  nation  tests  on  (+45/902)s  specimens 
were  conducted  to  study  the  interaction  among  matrix  degradation,  edge  crack¬ 
ing/ply  delamination  and  intraply  cracking.  This  was  done  jointly  with  the 
Mechanics  and  Materials  Center  at  Texas  A&M  University  where  they  have  a  simi¬ 
lar  interest  on  the  subject  under  an  Air  Force  contract  F49620-82-C-0057.  The 
Phase  IV  program  provided  Texas  A&M  University  ten  residual  strength  specimens 
and  thirty-six  delaimination  specimens  for  testing.  In  return,  Texas  A&M 
University  provided  the  Phase  IV  program  results  of  their  findings. 


During  testing  at  A&M  University,  strain  measurement  was  made  by  using  an 
extensometer  called  an  "Elastic  Strain  Ratio  Translater"  (ESRT)  gage. 

This  ESRT  extensometer  was  developed  by  B.  C.  Harbert  of  A&M  University  and 
its  description  is  shown  in  Figure  23.  Vought  used  the  extensometer  with  two 
LVDTs  as  described  in  Figure  9. 


Tko  f  or+  iw'A  f  Aw»  +Via  r  tw*Anri4*F»  4- nr  fr  r  a  Um.in  i  n  C-J  hiima  O  A 

•  •  i  —  vv  j  v  pi  vvuuui  v  i  \j  i  v»i  i  v.  i  uj  •  uuu  i  on  uc.  j  oa  ■  a  ouonn  *  i  ■  i  i  vju  i  v;  c-“t  • 

Although  Vought  and  A&M  used  slightly  different  procedures,  it  involved  essen¬ 
tially  the  same  test  steps:  static  preload,  fixed  duration  of  fatigue  cycling 


and  gradually  increased  static  loading  to  failure.  Tests  were  conducted  at 
75°F/50%  R.H.  and  132°F/50%  R.H.  environments  and  two  fatigue  stress 


levels,  55%  of  F^,  and  60%  of  .  Test  results  obtained  by  Vought  are 
summarized  in  Tables  14  and  15.  Fatigue  cycling  was  shown  to  have  degrada- 
tional  effects  on  material  stiffness  and  strength.  It  does  not  seem  to  affect 
the  failure  strain.  Edge  cracking  was  monitored  during  the  residual  strength 
test  by  using  thin  replicating  tape.  The  typical  positive  prints  from  the 
replication  technique  are  shown  in  Figure  25.  Transverse  cracks  across  the 
four  90°  layers  and  delamination  along  the  interface  betwen  the  90°  layer 
and  45°  layer  are  the  two  common  crack  formations  along  the  edges.  The  rep¬ 
lication  technique  was  used  to  monitor  the  cracking  formation  and  results  are 
shown  in  Tables  16  and  17.  At  low  fatigue  stress  levels  (55%  of  F^u),  trans¬ 
verse  cracking  occurs  earlier  than  delamination.  At  a  higher  stress  level  (60% 
of  Ff-ij).  the  delamination  cracks  occur  first. 
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FIGURE  24  RESIDUAL  STRENGTH  TEST  PROCEDURE 
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DEVELOPMENT  OF  CROCKS  FOR  (+45/90-,)  SPECIMENS  AT  75  F/50%  RH  ENVIRONMENT 


DEVELOPMENT  OF  CRACKS  FOR  (+45/90.,),  SPECIMENS 


Residual  strength  tests  were  also  conducted  by  A&M  University  for  regular 
(H5/902)s  specimens  and  delamination  specimens  described  in  Section  3.0. 

Their  test  results  are  shown  in  Table  18.  The  typical  data  plots  are  shown  in 
Figure  26  based  on  the  test  procedure  AM-A  type  test  of  Figure  24.  The  hys¬ 
teresis  loops  in  the  stress-strain  curves  are  not  obvious  after  the  fatigue 
cycling  and  gradually  open  up  under  gradually  increasing  static  loads.  It  was 
found  from  the  results  in  Table  18  that  center  defects  had  a  significant 
affect  on  the  specimens'  modulus  and  strength. 
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+  REFER  FIGURE  24  FOR  TEST  TYPE 
>  RELATIVE  HUMIDITY  NOT  CONTROLLED 


12??4 


tti*?  .  00382 


31' :  CYCLING 


.00009 


•00083  .0015*  .0023  .00304  ,00j?8 


STMIN 


Cb)  RAMP  TO  55X  AFTER  FATIGUE  CYCLING 


**  >  ’  i  ” 

....  i 

l  Kf  j  i  I 


K 


(c)  RAI 


5  .0051? 


Reproduced  From 
Best  Available  Copy 


y  -y 

,4*  ^ 


sV 


.0000?  nnii  1  ~1 

Wn  .00014  .0051?  .0042  .00523  .00836 


strain 


(e)  RAMP  TO  80S!  AFTER  FATIGUE  CYCLING 


\C. 


(f)  RAMP 


FIGURE  26  TYPICAL  STRESS-STRAIN  CURVES  AT  VARIOUS  STAGES  OF  THE  RESIDUAL  STRENGTH  TEST 
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STAGES  OF  THE  RESIDUAL  STRENGTH  TEST 


7.0  ANALYSIS  OF  DAMAGE  GROWTH  AND  FAILURE 


In  the  Phase  III  report,  a  theoretical  model  based  on  matrix  degradation 
was  proposed  for  estimating  the  uniaxial  laminate  failure  strain  in  static  and 
fatigue  loading  and  for  predicting  fatigue  life.  Specifically,  it  was  shown 
that  by  using  a  matrix  modulus  of  approximately  10%  of  its  initial  (undamaged 
state)  value,  theoretically  predicted  strains  were  close  to  those  at  failure 
in  static  and  fatigue  tests.  Although  many  of  the  failed  coupons  exhibited 
considerable  delamination,  it  was  suggested  that  any  significant  amount  of 
delamination  occurred  primarily  in  the  latter  stages  of  failure  and  was  trig¬ 
gered  by  gross  matrix  damage  within  plys. 

These  ideas  have  been  further  investigated  as  part  of  the  Phase  IV 
effort.  In  this  section  we  consider  in  more  detail  the  separate  effects  of 
matrix  damage  and  delamination,  and  indicate  possible  ways  in  which  they 
interact  to  produce  overall  failure.  Since  writing  the  Phase  III  report,  some 
related  studies  have  been  published  by  others10'13  which  are  very  informa¬ 
tive  on  the  nature  of  intraply  damage  and  delamination  in  laminates  similar  to 
ours.  This  information,  as  well  as  some  recent  theoretical  work,H,/  will 
be  taken  into  account  in  this  section. 

7.1  DAMAGE  THEORY  FOR  LINEAR  ELASTIC  COMPOSITES 

According  to  the  theory  of  linear  elastic  fracture  mechanics,  growth  of  a 
crack  initiates  when  the  so-called  energy  release  rate  G  exceeds  a  critical 
value,  Gc .  This  critical  value  is  material -dependent  and,  in  general, 
depends  on  environmental  factors  such  as  temperature  and  moisture  content  in 
the  case  of  plastics.  Thus,  an  important  aspect  of  damage  growth  analysis  is 
the  prediction  of  energy  release  rates  when  cracks  represent  at  least  part  of 
the  damage.  As  there  is  considerable  evidence  that  any  significant  damage  in 
advanced  polymeric  composites  usually  consists  of  intraply  and  delamination 
cracks,10"13  we  shall  discuss  the  prediction  of  G  for  both  types  of  crack¬ 
ing.  Extension  of  the  results  to  viscoelastic  behavior  is  then  made  using 
existing  theory. 
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Consider  a  tensile  specimen  in  which  the  applied  axial  load  is  F  and  the 
associated  axial  displacement  between  the  points  of  load  application  is  U. 
Assuming  linear  elastic  behavior  for  a  fixed  state  of  damage, 

F  =  kU  (4) 


where  k  is  the  stiffness.  In  general,  k  is  a  function  of  damage,  which  shall 
be  represented  by  a  set  of  N  crack  surface  areas  Ai  (i=l, — N).  Neglecting 
residual  strain  energy,  the  total  strain  energy  in  the  coupon  is 

W-,  =  1/2  kU2  (5) 

where  k  =  k ( ) -  The  familiar  equation  for  energy  release  rate  is  obtained 
by  forming  the  derivative  of  Wy  with  respect  to  A^ , 


1  3k  „2 

'  7  3A~  U 


(6) 


where  k  -  k  ( Ay )  and  i,  j  -1  ,2 , — N.  Physically,  G.jdA^  (i  not  summed)  is 

the  work  done  by  the  elastic  material  on  the  failing  materia1  at  the  ith  crack 

tip  where  the  area  of  this  crack  increases  by  an  infinitesimal  amount  dA^. 

Inasmuch  as  U  is  constant  when  this  change  is  evaluated,  the  externally 

applied  forces  do  not  work,  and  the  entire  change  in  strain  energy  dWy  is 

con. jrted  to  mechanical  work  at  the  advancing  crack  tip.  Strictly  speaking, 

the  notation  which  characterizes  the  shape  and  orientation  of  each  existing 

14 

crack  and  newly  formed  area  should  be  used  .  However,  for  present  purposes 
it  will  be  sufficient  to  use  a  simplified  notation  in  which  only  the  magni¬ 
tudes  of  the  areas  A^  are  explicitly  identified. 

Let  G^..  be  the  critical  value  of  ,  i.e.,  is  the  "energy"  (more 
precisely,  the  "mechanical  work")  per  unit  area  required  to  separate  material 
points  at  the  crack  tip  associated  with  the  new  area  dA^ .  Whether  or  not 
the  new  bu^face  dA^  is  actually  formed  depends  on  whether  or  not  the  avai 1  - 
ahl  e  energy  exceeds  the  required  energy  Gc  ^ ;  i.e.,  if  <GC1-  then 
dA^  *  0  and  if  (i-  :  Gci  then  dA^  >  0. 
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By  defining  compliance,  C  =  i/k,  and  using  the  relation  U  =  CF,  equation 
(6)  may  be  rewritten  in  a  similar  form  in  terms  of  applied  load. 


6 


(7) 


Equations  (6)  and  (7)  are  standard  forms  for  the  energy  release  rate  of  speci¬ 
mens  and  structures  when  a  single  applied  load  and  displacement  are  consid¬ 
ered.  In  predicting  for  specimens  under  original  stress,  it  is  often 
helpful  to  rewrite  these  equations  in  terms  of  effective  axial  modulus,  E, 
overall  or  average  stress  and  strain  e  .  Thus, 

with  o  =  Ee,  a  =  F/Acs  and  e  =  —  (where  =  cross-sectional  area  and  L  = 
length  of  coupon  prior  to  application  of  F),  we  obtain 


k  =  EAcs/L,  C  =  LU/AC$ 


where  D  =  E  ^  is  the  uniaxial  compliance.  Thus, 


V  3E  2 
=  ~  I  *7T  r 

i 


or,  equivalently, 


i  2 


V  30  2 

0 


in  which  V  =■  is  the  volume  of  the  sped 


men. 


(8) 


(9) 


(10) 


In  studies  of  the  graph i te/ epoxy  composite  material  T300/5208,  which  is 

11  l  ? 

very  similar  to  the  composite  used  here,  O'Brien  ’  found  that  the  modulus 
E  of  the  laminates  decreased  linearly  with  the  amount  of  delamination  area 
(as  defined  by  the  area  projected  onto  the  plane  of  the  laminate  through 
x-ray  measurements).  Thus,  using  O'Brien's  notation  for  moduli. 


1  '  <F*-W  A 


dT 


'LAM 


(11  ) 
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where  A^y  =  2bL  is  the  projected  delaminated  area  of  a  fully  delaminated 
sample  of  length  L.  Noting  that  Ad  =  0  without  del  ami  nation,  and  Ad  = 
A^y  with  full  delamination,  equation  (11)  implies  that 


ELAM  “  sPecimen  modulus  without  delamination 

E*  =  modulus  of  a  fully  delaminated  specimen 

As  an  illustration  of  Ad  and  A^y,  consider  the  laminate  in  Figure  27, 
which  shows  the  end  view  of  a  specimen  with  edge  del  ami  nations  cf  depth  A.  If 
the  depth  were  constant  along  the  specimen  length  L  (x-direction) ,  then  A^  = 
2AL.  O'Brien  assumed  a  pair  of  edge  delaminations  existed  above  and  below  the 
center  plane  of  symmetry  in  analytical  predictions  of  modulus;  therefore,  four 
crack  tips  (rather  than  the  two  in  Figure  27)  were  assumed  to  be  involved  in 
the  del  ami  nation  process.  However,  Ad  and  AdT  are  the  same  for  the  cases 
of  two  and  four  cracks  in  view  of  their  definition. 

Recall  that  the  energy  release  rate  G_.  is  the  mechanical  work  available 
at  a  crack  tip  per  unit  of  new  area  at  this  one  tip:  thus,  only  one  crack  tip 
is  advanced  when  calculating  the  derivative  E/  Ai  in  Equation  (9).  with 
this  point  in  mind,  let  us  first  consider  the  two-crack  case  in  Figure  27,  in 
which  we  let  i  =1  and  i  =  2  represent  the  two  cracks,  and  E*^  the  specimen 
modulus  with  full  delamination.  Using  Ad  =  A-j  +  ^  ?nd  Ady  =  2bL  in 
Equation  (11)  for  E,  we  then  obtain  from  Equation  (9)  for  i  =  1, 


G 


i 


(E 


LAM 


-  e2  ) 


(12) 


In  the  case  of  fot.'*"  cracks  where  the  cracks  extending  from  the  same  edoe  have 

the  same  depth,  one  has  A^  =  ^  +  A^  +  A3  +  A^)/V  and  A^y  :  2bl  ; 

thus 


G ,  = 


11. 


(E 


LAM 


* 

t4  ) 


^  !l|u,,-ius  when  the  specimen  is  fully  del 


ami na ted. 


(13) 
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The  value  of  energy  release  rate  is  the  same  for  each  of  the  two  cracks 
(Equation  (12))  or  each  of  the  four  cracks  Equations  (13).  Furthermore,  these 
two  equations  are  really  not  restricted  to  equal  or  uniform  crack  depths. 
Namely,  each  crack  could  have  a  different  depth  and  it  could  vary  along  the 
length;  indeed,  lengthwise  variation  was  observed.^  Although  the  formula 
for  G.  is  independent  of  depth  (because  E  is  linear  in  delamination  area), 
it  does  depend  on  the  number  of  delamination  cracks  (through  factor  1/2  or  1/4 
and  the  value  of  E*2  or  E*^).  O'Brien's  predictions  were  based  on  Equa¬ 
tion  (12),  but  with  E*^  used  in  place  of  E*2.  This  resulted  in  the  best 
agreement  for  measured  stiffness  changes  and  the  critical  value  of  G^.  The 
need  for  this  modification  of  the  basic  formula  for  G-j  is  probably  due  to 
the  presence  of  many  cracks  within  each  ply,  which  were  not  accounted  for  in 
predicting  and  E*,  and  their  participation  in  forming  the  delamination 
crack  tips;  these  intralaminar  cracks  are  shown  and  discussed  in  detail  in 
References  10  and  11.  O'Brien  found  that  essentially  only  one  delamination 
grew  from  each  edge,  but  it  skipped  from  one  side  of  the  center  plane  to  the 
other  through  transverse  cracks. 

The  effect  of  intralaminar  cracks  on  both  E.  ...  and  E*  will  be  considered 

LHI'I 

here.  First,  however,  we  shall  make  some  comparisons  using  Equations  (12)  and 
(13)  which  provide  evidence  for  the  fact  that  static  and  fatigue  failure 
behavior  reported  on  the  program  cannot  be  explained  without  accounting  for 
th  em . 

The  modulus  ELAM  without  damage  for  the  [+45/90,,]s  laminate  has  been 
predicted  elsewhere  in  this  report  (Figure  11);  the  largest  and  smallest  val¬ 
ues,  based  on  the  extremes  of  constituent  properties  for  all  environments 
studied,  are  3.38  and  3.17  x  106  psi.  The  modulus  with  full  delamination  E* 

will  be  predicted  using  classical  lamination  theory  for  several  cases;  Jones' 

1 5 

notation  is  employed.  For  a  laminate  in  terms  of  the  coordinates  in  Fig¬ 
ure  27,  the  forces  N  =  (Nx,  Ny,  Nxy),  moments  M  =  (Mx,  M  ,  Mxy), 
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1 


middle  surface  strains  £°  =  (^x°>  £y°.  )  and  curvatures  <  -  (kx,  k^) 

are  related  according  to 


(14) 


and,  in  the  inverse  representation, 

§if]  r  n  i 

b”  d~ j  |tJ 


(15) 


The  6  x  6  symmetric  matrices  in  Equations  (14)  and  (15)  are  obviously  inverses 
of  one  another,  and  are  each  comprised  of  three  3x3  symmetric  matrices  (A, 

B,  D)  and  (A1,  B',  [)'). 

The  modulus  for  uniaxial  loading  without  delamination,  ELAfJ|,  is 
elam  (h  A<  i  r 

where  A'^  is  the  first  element  of  the  compliance  matrix  A'.  The  modulus 
with  complete  delamination,  E*,  in  which  all  of  the  n  separated  laminates 

(sub laminates)  are  under  the  same  axial  strain,  is 

o 

=  K  C  Eihi  U7) 

i=l 
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By  definition, 


E1  *  Wx 


08) 


which  is  the  axial  modulus  of  the  i th_  sublaminate,  h^  is  its  thickness,  and 
NX1-  is  the  axial  force  on  it. 


It  is  not  obvious  what  constraints  or  edge  loadings  should  be  used  in  pre¬ 
dicting  N  O'brien^  assumed  all  three  curvatures  k  and  edge  forces 
Ny.j  and  NXyj  are  zero;  this  may  be  an  acceptable  application  for  some 
cases  in  view  of  the  constraint  imposed  on  out-of-plane  deformations  by  the 
uncracked  center  section  along  the  delamination  crack  tips  and  the  alternating 
plane  of  delamination  along  the  sample  length  observed.  For  zero  curvatures 

S"11  Nyi  =  Nxyi  =  °> 


Ei  * 


-1 


(19) 


where  A"^  is  the  first  element  of  the  inverse  of  the  submatrix  for 
the  itlh  sublaminate.  For  relatively  deep  delaminations  without  any  or  many 
changes  in  the  plane  of  delamination,  significant  curvatures  of  unbalanced 
and/or  unsymmetric  sublaminates  may  occur.  As  an  extreme  case,  where  each 
sublaminate  is  under  no  edge  moments  and  forces  other  than  Nx..,  the  axial 
modulus  is 


E. 


<hiA'ni  > 


(20) 


where  A'^-  is  the  first  element  of  the  matraix  A'  for  the  ith  sublaminate. 


The  moduli  (19)  and  (20)  enable  estimates  of  the  energy  release  rate  to  be 
made.  In  order  to  obtain  more  accurate  predictions,  one  should  analyze  the 
partially  delaminated  test  specimen  using  plate  theory  or  a  more  general 
approach,  allowing  for  spacially  varying  strains  and  curvatures.  However, 
this  much  more  complicated  problem  will  not  be  addressed  here.  Instead,  we 
shall  consider  one  more  case,  which  is  intermediate  to  the  zero  curvature 
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result.  Equation  (19),  and  unconstrained  sublaminates,  Equation  (20).  It  is 
motivated  by  the  problem  In  which  delamination  occurs  along  the  center  plane 
of  the  C+45/902]s  laminate.  Assuming  delaminations  of  equal  depth  extend 
from  each  edge,  and  that  the  depth  is  constant  along  the  sample  length,  the 
zero  curvature  case.  Equation  (19)  yields  E-  =  E^;  thus,  the  energy 
release  rate  is  predicted  as  zero.  In  practice,  an  opening  mode  of  crack  dis¬ 
placement  occurs,  in  which  the  curvature  serves  to  define  the  amount  of 
opening  displacement.  For  a  relatively  long  specimen,  a  reasonable  assumption 
is  that  the  out-of-plane  displacements  do  not  vary  along  the  length,  (K^  = 

Kx  =  0),  and  the  free  edge  forces  and  moments  are  zero  (Ny  =  =  My 

=  =  0).  Inasmuch  as  Ex  is  specified  to  be  the  same  for  all  sublami¬ 

nates,  one  too  many  conditions  have  been  specified;  we  shall  drop  the  require¬ 
ment  =  0  in  one  problem,  and  drop  Nxy  =  0  in  a  second  problem. 


Moduli  for  the  [+45/90  9  ]  c  laminate:  The  center-plane  delamination 
problem  will  be  discussed  first.  For  the  constituent  properties  which  result 
in  the  largest  modulus,  ELW  =  3.38  x  106  psi ,  we  find 

E*  =  2.42  x  106  psi  if  Ny  =  Nxy  =  My  =  Kx  =  i^y  =  0  (21 ) 

and  only  a  slightly  different  value 

E*  =  2.48  x  106  psi  if  Ny  =  Mxy  =  My  =  Ky  =  l^y  =  0  (22). 


Equation  (21)  provides  a  higher  energy  release  rate  than  (22),  and  will  be 
used  in  delamination  predictions.  For  the  fully  unconstrained  case.  Equation 
(20), 

E*  =  2.11  x  106  psi  (23) 

Referring  now  to  the  problem  illustrated  in  Figure  27,  in  which  one  pair 
of  cracks  exists  between  the  -45  and  90  plys,  we  find  for  the  properties  cor¬ 
responding  to  ELflM  =  3.38  x  106  psi , 

E*  =  2.80  x  10®  psi  (based  on  Equation  (19))  (24) 
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and 


(25) 


E*  =  2.02  x  106  psi  (based  on  Equation  (20) 

Next,  in  addition  to  the  cracks  shown  in  Figure  27,  assume  another  pair  of 
cracks  exist  between  the  90  and  -45  plys  below  the  sample  middle  surface.  Then 

E*  =  2.03  x  10^  psi  (based  on  Equation  (19)  (2b) 

and 

E*  =  1.85  x  106  psi  (based  on  Equation  (20)  (27) 


Critical  axial  strain  for  delamination  of  the  [t45/90?]c _ 1  ami nat e: 

Equations  (12)  and  (13)  may  be  used  to  predict  the  strain  at  which  delamina- 
tion  occurs,  e  c,  given  the  critical  value  ec  : 


,  2Gc  ,0.5 


(2  cracks) 


(28) 


4G. 


0.5 


(4  cracks) 


(29) 


where,  by  definition 


AE  = 


elam 


E* 


(30) 


For  the  opening  mode  of  delamination  in  the  AS 1 / 3501 - 6  composite  in  moder- 

1  ft 

ate  environments,  Gc  M  lb/in  as  reported  by  Wilkens.  Mixed  mode  delami¬ 
nation  G  is  higher,  and  a  value  of  G  2  ]  b/ i n  will  be  used  here.  (The 
C  Q 

measurements  of  G  c  for  mixed-mode  cracking  were  made  on  unidirectional  0 
fiber  angle  laminates,  and  it  may  be  different  for  the  laminate  under  study 
here.)  Strains  predicted  from  Equations  (28)  and  (29)  are  listed  below,  in 
which  the  number  in  parenthesis  corresponds  to  the  particular  E*  which  was 
used  in  At;  also,  Elm  =  3.38  x  106  psi  and  h  =  .044  inches. 
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NUMBER  AND  LOCATION 


CASE 

CRITICAL  STRAINS  (%) 

E*  EQUATION 

Gc  (Ib/in) 

OF  CRACKS 

1 

.69 

(21) 

1 

2,  90/90 

2 

.60 

(23) 

1 

2,  90/90 

3 

.89 

(24) 

1 

2,  -45/90 

4 

1.26 

(24) 

2 

2,  -45/90 

5 

.58 

(25) 

1 

2,  -45/90 

6 

.82 

(25) 

2 

2,  -45/90 

7 

.82 

(26) 

1 

4,  -45/90 

8 

1.16 

(26) 

2 

4,  -45/90 

9 

.77 

(27) 

1 

4,  -45/90 

10 

1.10 

(27) 

2 

4,  -45/90 

When  the  low  values  of  constituent  moduli  are  used,  corresponding  to 
fS  ^ 

E.  =  3.17  x  10  psi,  the  critical  strains  are  predicted  to  be  less  than 
LAM 

those  in  the  table,  but  only  by  a  small  amount;  the  ratio  of  these  strains  to 
those  in  the  table  is  approximately  0.96. 

It  is  instructive  to  compare  the  critical  strains  to  those  at  which  stiff¬ 
ness  changes  have  been  reported  in  the  experimental  studies  for  static  and 
fatigue  loading.  In  general,  all  of  the  predicted  critical  values  are  above 
those  at  which  specimens  begin  to  exhibit  noticeable  softening  due  to  damage. 
Only  for  cases  2  and  5  do  the  strains  approach  the  experimental  values;  but  it 
is  believed  these  predictions  are  unrealistically  low  as  the  sublaminates  were 
assumed  to  be  completely  unconstrained.  Also,  except  for  samples  with  brass 
inserts  at  the  middle  plane,  no  sigificant  delaminations  at  the  middle  surface 
were  observed  in  failed  coupons.  In  view  of  these  observations,  we  estimate 
from  the  theory  that  £  ^  .7%  for  the  samples  with  centered  brass,  and 

V  ' 

£  <1.8%  for  the  other  samples;  most  delaminations  appeared  to  occur  between 
the  -45  and  90°  plys,  regardless  of  whether  or  not  brass  separators  were 
used  between  these  plys.  It  should  also  be  noted  that  the  energy  release  rate 
during  fatigue  delamination  of  0°  laminates  was  close  to  the  static  critical 
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value  Gc,  and  therefore  we  have  no  basis  at  this  time  for  using  in  Equa¬ 
tions  (28)  and  (29)  an  energy  release  rate  which  is  significantly  less  than  1 
lb/in. 

Intraply  Damage  and  Viscoelastic  Effects:  On  the  basis  of  the  above  esti¬ 
mates  and  the  considerations  in  Phase  III  concerning  the  effect  of  a  reduced 
matrix  modulus,  we  believe  intraply  microcracking  has  to  be  taken  into  account 
besides  delamination  in  order  to  predict  the  deformation  and  failure  behavior 
of  laminates;  indeed,  intraply  damage  may  account  for  the  largest  amount  of 
softening  prior  to  failure  in  our  specimens,  considering  how  large  the  pre¬ 
dicted  critical  strains  without  damage  are.  Measurements  of  damage  using 
X-rays  would  be  helpful  to  confirm  the  point.  O'Brien11  predicted  that 
delamination  caused  more  softening  than  intraply  cracking  in  the  laminates  he 
studied,  and  this  was  supported  by  direct  measurement  of  the  amount  of  delami¬ 
nation. 


The  intralaminar  cracking,  as  defined  by  parameters  ,  reduces  E,  Equa¬ 
tion  (11),  in  which  both  E,  and  E*  are  reduced  by  this  cracking. 


14 


bcnapery' '  developed  a  relatively  general  theory  for  elastic  and  viscoelas¬ 
tic  materials  (allowing  for  a  distribution  of  and  values)  which 
could  be  used  to  account  for  the  effect  of  strain  history  and  viscoelastic 
properties  on  the  damage  and,  in  turn,  its  effect  on  moduli.  However,  explic¬ 
it  realistic  predictions  of  damage  are  very  involved,  and  will  not  be  done 
here.  This  complexity  exists  whether  or  not  viscoelastic  effects  are  consid¬ 
ered.  For  example,  the  stiffness  E*  depends  on  not  only  the  damage  which 
occurs  in  the  sublaini nates  after  they  are  formed,  but  also  on  the  damage  prior 
to  delamination  and  that  produced  by  the  high  stresses  around  a  delamination 
crack  tip  v/hile  it  is  moving  and  producing  sublaminates.  In  the  center  part 
of  a  sample  which  is  not  delaminated,  another  complication  exists  because  of 
the  breakdown  of  classical  lamination  theory  (CIT).  Arc-nburg^  has  shown 
that  when  the  effective  matrix  modulus  Em  is  reduced  from  its  initial  value 
of  0.6  x  106  psi,  interlaminar  shear  deformations  may  become  important.  For 
the  [+_45/90(i]s  specimens  used  in  this  program,  energy  release  rates  and 
stiffness  predicted  by  CLT  are  iri  considerable  error  for  damage  states  at  the 
failure  strains.  This  behavior  seems  to  imply  one  should  account  for  strain 
nonuniformity  in  the  undelaminated  part  of  the  laminate  in  order  to  predict 
damage.  Complex  temperature  dependence  of  Gc  for  mixed-mode  cracking16 
and  residual  stresses  further  complicate  the  problem. 
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8-0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  effect  of  temperature  and  moisture  on  the  deformation  and  failure 
behavior  of  AS/3501-6  uniaxial  lamina  and  laminates  is  very  complex.  The  com¬ 
monly  employed  assumption  of  thermorheologically  simple  behavior  for  polymeric 
materials  (in  which  changes  in  temperature  and  moisture  simply  alter  the 
time-scale  for  response)  is  not  valid.  However,  the  Phase  I  studies  indicated 
that  creep  and  recovery  strains  of  matrix-dominated  specimens  with  a  constant 
damage  state  can  be  characterized  using  a  slightly  generalized  representation; 
namely j  besides  the  environmental  effect  in  the  response  time-scale,  the  envi¬ 
ronment  affects  the  initial  creep  compliance.  When  the  damage  state  changes 
continuously,  as  in  the  fatigue  tests,  it  has  not  been  possible  to  account  for 
temperature  and  moisture  effects  in  any  simple  way.  The  S-N  curves  for  a 
given  lamina  or  laminate  can  be  approximated  by  a  straight  line  using  log-log 
axes,  and  in  some  cases  (e.g.,  the  [+45/902]s  laminate)  the  slope  is  the 
same  for  all  environments  studied.  However,  in  general,  the  slope  (i.e.,  the 
power  law  exponent)  is  different  for  different  laminates.  Also,  in  some  cases 
an  increase  in  temperature  and/or  moisture  reduces  the  fatigue  life,  where  in 
other  cases  it  increases  the  fatigue  life.  It  is  believed  this  complex  fail¬ 
ure  behavior  is  due  at  least  in  part  to  the  delamination  fracture  toughness 
for  the  opening-mode  increasing  with  increasing  temperature,  and  for  a  mixed 
mode  (with  a  large  shearing  component)  exhibiting  the  opposite  behavior.16 

Intralaminar  cracks  appear  to  have  a  significant  effect  on  the  initiation 
and  growth  of  delamination  cracks.  Theoretical  estimates  of  energy  release 
rates  available  to  drive  delaminat.ions  give  values  which  are  considerably  below 
those  required  for  their  propagation  unless  interlaminar  cracks  are  taken  into 
account.  It  is  thus  believed  that  both  types  of  damage  have  to  be  considered 
to  predict  deformation  behavior  of  laminates  with  damage  and  to  predict 
fatigue  lifetimes.  An  approach  using  results  from  viscoelastic  fracture 
mechanics  theory  is  discussed.  But.,  in  order  to  develop  and  verify  a  realis¬ 
tic  theoretical  model,  it  appears  necessary  to  make  extensive  use  of  NDI  meth¬ 
ods  at  various  stages  of  specimen  damage  in  order  to  identify  the  actual  dam¬ 
age  modes,  A  study  of  new  specimen  geometries  and  fiber  layups,  taking  into 
account  data  already  generated,  should  be  made  in  order  to  separate  as  much  as 
possible  distinct  effects  of  intralaminar  cracking,  del  ami  nation,  and  residual 
stresses. 
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